Double Extensive Air Shower Induced by Ultra-High Energy Cosmic Tau-Neutrino 
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We investigate the possibility of detecting ultra-high energy cosmic tau-neutrinos by means of 
a process involving a double extensive air shower, the so-called Double-Bang Phenomenon. In 
this process a primary tau-neutrino interacts with an atmospheric quark creating a first hadronic 
shower and a tau-lepton, which subsequently decays creating a second cascade. The number of these 
events strongly depends on the flux of tau-neutrinos arriving at the Earth's atmosphere and can be 
used to test some theoretical models related to the production of ultra-high energy tau-neutrinos. 
We estimate the potential of the fluorescence detector of the Pierre Auger Observatory to observe 
Double-Bang events. We conclude that for tau-neutrinos with energies ranging from O(O.l) EeV 
to O(10) EeV the number of detected events vary from hundreds in a year to only few events in 
hundreds of years. 

PACS numbers: 13.15.+g, 96.40.Pq 



I. INTRODUCTION 

It is believed that ultra-high energy cosmic neutrinos 
may play an important role to explain the origin of cos- 
mic rays with energies beyond the GZK limit of few times 
10 19 eV pl@, once that neutrinos hardly interact with 
cosmic microwave background or intergalactic magnetic 
fields, keeping therefore its original energy and direction 
of propagation. Even if they have masses or magnetic 
moments, or travel distances of the order of the visible 
universe, those characteristics do not change very much. 
Possible sources of these ultra-high energy neutrinos, like 
Active Galactic Nuclei and Gamma Ray Bursts, are typ- 
ically located at thousands of Mpc 0, Considering 
that neutrinos come from pions produced via the pro- 
cess 7+p^A r + 7r0], that there is an additional v e 
flux due to escaping neutrons and that about 10% of the 
neutrino flux is due to proton-proton (pp) interactions, 
the proportionality of different neutrino flavors result: 
v e : : v T = 0.6 : 1.0 :< 0.01 Nevertheless, obser- 
vations of solar and atmospheric neutrinos present 
compelling evidence of neutrino flavor oscillations. Such 
oscillations have been independently confirmed by terres- 
trial experiments. KamLAND |8( observed D e disappear- 
ance confirming (assuming CPT invariance) what has 
been seen in solar neutrino detections and K2K 
served v^jv^ conversion compatible with what has been 
detected in atmospheric neutrino observations. 

In order to understand these experimental results 
by means of neutrino oscillations, two scales of mass 
squared differences and large mixing angles have to be 
invoked. For solar and KamLAND observations, AttIq ~ 
7 x 10~ 5 eV 2 and sin 2 20 Q ~ 0.8. And for atmo- 
spheric neutrino and K2K, |Arn 2 tm | ~3x 10 -3 eV 2 and 
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sin 20 a t m ~ 1- Moreover LSND experiment |TT| may 
have observed 9^ — ► v e transition which can be also ex- 
plained by neutrino oscillations with a large mass scale, 
\&m 2 LSND \ ~ (0.5 - 2.0) eV 2 . Such results will soon be 
checked by MiniBooNE 12]. These scales require four 
neutrino oscillation framework (or three, if LSND re- 
sults will not be confirmed by MiniBooNE experiment) 
which imply, for ultra-high energies of the order 1 EeV 
or higher, oscillation lengths much smaller than typical 
distances from the sources of ultra-high energy neutrinos. 
Consequently when neutrino flavor oscillations are taken 
into consideration the flavor proportion will be modified 
to v e : j/ju : v T ~ 1 : 1 : 1. Therefore one expects a con- 
siderable number of tau-neutrinos arriving at the Earth. 

In this paper we investigate the possibility of detect- 
ing ultra-high energy cosmic tau-neutrinos by means of a 
process in which a double Extensive Air Shower (EAS) is 
identified, the so-called Double-Bang (DB) Phenomenon. 
In that kind of event a tau-neutrino interact with a quark 
via charged current creating one cascade of hadronic par- 
ticles and a lepton tau which subsequently decays pro- 
ducing a second cascade. DB Phenomenon was first pro- 
posed for detectors where the neutrino energy should be 
around 1 PeV [f| . It does not happen with neutrinos dif- 
ferent from tau. The electron generated by an electron- 
neutrino immediately interacts after being created and 
the muon generated by a muon-neutrino, on the other 
hand, travel a much longer distance than the size of the 
detector before interacting or even decaying. So we do 
not have DB events from them. 

In order to identify a DB Phenomenon we have to look 
for two Extensive Air Showers (EAS) in the same direc- 
tion of propagation inside the field of view (f.o.v.) of the 
detector, i.e., in the physical space around the detector 
in which an event can be triggered. 

In the Pierre Auger Observatory a hybrid detec- 
tion technique will be used to make a detailed study 
of cosmic rays at energies mainly around 10 EeV and 
beyond. The two techniques consist of an array of de- 
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tectors spread on the ground (the ground array detec- 
tor also called surface detector) and an optical detec- 
tor used to probe longitudinal development of EAS by 
recording the fluorescence light emitted by the excited 
nitrogen molecules of the Earth's atmosphere. It has 
been shown that the Auger Observatory can detect at- 
mospheric near-horizontal air showers generated by neu- 
trinos with the surface detector [lj, [l5|, Il6l Il7l Il8l Il9| . 
We concentrate here on a different approach to study the 
possibility of detecting events induced by ultra-high en- 
ergy tau-neutrinos, the DB events, using the fluorescence 
detector of that observatory. 

We conclude that the features of the Pierre Auger Flu- 
orescence Detector favor the observation of DB events 
with tau-neutrino energies varying from 0(0.1 EeV) to 
0(10 EeV), despite the low efficiency of the fluorescence 
detector at energies smaller than 1 EeV. We estimate 
the number of DB events observed in the fluorescence 
detector varying from hundreds in a year to few events 
in hundreds of years depending mainly on the primary 
tau-neutrino flux. 

This paper is organized in the following way: Sect ion ITT1 
has a brief introduction to the DB Phenomenon. Sec- 
tion [^21 shows some results of the DB events simulated 
and Section lll Bl describe how we calculate the number of 
events in the Pierre Auger Fluorescence Detector. Sec- 
tion |n3]has the number of events calculated for different 
models of ultra-high energy neutrino flux and discuss how 
could we take some physical information from that. The 
conclusions are in Section IIVI in which we discuss also 
the background events. 



II. THE ULTRA-HIGH ENERGY 
DOUBLE-BANG AND THE AUGER 
OBSERVATORY 

Studying the characteristics of the fluorescence detec- 
tor, such as its efficiency and f.o.v. and the character- 
istics of the DB events generated by ultra-high energy 
tau-neutrinos, one can estimate the rate of that kind of 
event expected in the Auger Observatory. 

Fig. n shows a schematic view of an ultra-high energy 
DB with the detector position and the time integrated 
development of the two showers, one created by the high- 
energy tau-neutrino interacting with a nucleon in the at- 
mosphere and the other created by the decay of the tau 
generated in the first interaction of the tau-neutrino. The 
f.o.v. of the fluorescence detector will be comprehended 
between angles near the horizontal (~ 2°) and 30°, and 
a radius r of approximately 30 km. The maximal height 
from where the DB can be triggered by the fluorescence 
detector is h and to is its projection in the DB propaga- 
tion axis. The zenith angle is represented by 9. 

We considered only showers moving away from the 
detector since, in the opposite large amount 

of Cerenkov light arrives together with the fluorescence 
light, spoiling a precise data analysis po|. 




FIG. 1: A schematic view of a Double-Bang and the f.o.v. of 
the Pierre Auger Fluorescence Detector. See text for details. 



The total amount of light emitted by the first cascade 
is related to the energy transfered to the quark at the mo- 
ment of the first tau-neutrino interaction, which we will 
define as E±. The neutrino energy E v is the sum of the 
tau energy E T and E\, i. e., E v = E\ + E T . For charged 
current interactions above 0.1 EeV, approximately 20% 
of the neutrino energy is transfered to the quark [2l| . 
The second cascade, resulting from the tau decay, car- 
ries an energy E2 of approximately 2/3 E T and is visible 
when the tau decay is hadronic, which happens with a 
branching ratio of around 64% [22| ■ 

Therefore, very roughly, we have (Ei) ~ l/bE v and 
(E 2 ) ~ 2/3 (E T ) w 8/15^ and the relation between E\ 
and E 2 is given by: E 2 /E 1 = ^E V /\E V m 2.67. The 
distance traveled by the tau before decaying in laboratory 
frame is L = jet, where 7 = E T /m T and t is the mean 
life time of the tau, that has an error of approximately 
0.4% 0. 



A. Simulations 

In order to infer the possibility of observing DB events 
in atmosphere and detecting them in the Pierre Auger 
Fluorescence Detector, we first simulate the longitudi- 
nal development of the showers. For simplicity, we nu- 
merically simulate DB events using protons as primary 
particles. We create two separated showers using the re- 
lations of energy (E 2 /E 1 ) and distance (L) between the 
two EAS of a DB in the way presented in the previous 
section. In fact, it has been argued that it is difficult 
to distinguish if a ultra-high energy EAS was created by 
a proton or a neutrino |U 0] > at least when the cross 
section is extrapolated from the Standard Model. An 
important difference is the probability of interaction in 
the atmosphere, that can be 10 -5 smaller for a vertical 
neutrino than for a proton. There are also some works 
that study the differences between the longitudinal de- 
velopment of EAS generated by p rotons, heavier nuclei 
and different neutrino flavors |24l |25| . In reference pl| 
they use CORSIKA+Herwig Monte Carlo simulations to 
have electron and muon-neutrinos as primary particles, 
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TABLE I: Parameters used in the numerical simulations with 
the CORSIKA program and primary neutrino energy of 0.5 
EeV. For a visualization of the geometrical parameters 0, h 
and u>, see Fig.Q The units are km except where it is pointed. 



e a 


1st. mt. b 


Lcos(9 c 


2nd. int. d 


h c 


h (g/cm 2 )' 


CO 9 


45° 


24.3 


14.3 


10 


11 


231 


15.6 


55° 


19.6 


11.6 


8 


9.5 


291 


16.6 


65° 


15.5 


8.5 


7 


7.7 


377 


18.2 


75° 


10.2 


5.2 


5 


5.5 


515 


21.3 



"Zenith angle measuring the shower incident axis inclination 
''Altitude of the first interaction in the atmosphere 
"Projection in the vertical axis of the distance the tau runs before 
it decays 
d Altitude where the tau decays 

"Maximum altitude from where the fluorescence light of the DB 
can be detected, for r = 30 km (see Fig. 
■'Same as h, but in units of g/cm 2 , also called atmospheric depth 
'Projection of h in the DB propagation axis 



TABLE II: Same as Table Q where the primary neutrino en- 
ergy is taken to be 1 EeV. 



e 


1st. int. 


L cos 6 


2nd. int. 


h 


h (g/cm 2 ) 




45° 


33.6 


28.6 


5 


11 


231 


15.6 


55° 


28.2 


23.2 


5 


9.5 


291 


16.6 


65° 


22.1 


17.1 


5 


7.7 


377 


18.2 


75° 


13.5 


10.5 


3 


5.5 


515 


21.3 



but not tau-neutrinos. So the simulations we have made 
to study DB events still are a good approximation. We 
discuss how to distinguish events that could masquerade 
tau-neutrino induced DB events in the conclusions. 

Two different approaches can be taken to evaluate 
the longitudinal development of a DB. At first we used 
Gaisser-Hillas parameterization [2(|. This method, nev- 
ertheless is not reliable to determine the depth of the 
shower maximum and the point of first interaction, spe- 
cially for arrival angles larger than 60°. At this point, the 
CORSIKA simulation seems to be a good approach |27| . 
Here we use its version 6.00. 

Fig. [3 and Fig. [3] present the main results of the sim- 
ulations with the CORSIKA program. They show the 
longitudinal development of the DB for different incident 
angles chosen among the simulated events. Table [I] and 
Table Ull show the parameters used in the simulations. In 
these tables, the points of first and second interactions 
are chosen fixed parameters. The distance between the 
first and second interactions, as we pointed out above, is 
in accordance with the tau mean decay length (L) in the 
laboratory frame. Because of the small cross section of 
the neutrino, the probability of interaction in the top of 
the atmosphere is about the same for any point. 



B. The number of events 

To calculate the possible number of events in the Pierre 
Auger Fluorescence Detector, first consider for simplicity 
one detector with a site seeing of 360°. Then we can write 
the equation: 

N events = J P det (E, r, 9, <p) dE dA dT dQ (1) 

where, 

P det (E,r,9,tp) = Pint (E,0)x P trig (E,r,6,<p) ( 2 ) 

is the probability of detection of a DB in the Auger 
Observatory, given by the product of the probability 
Pnt{E,9) of the tau-neutrino to interact in the atmo- 
sphere and the probability P tr i g (E,r,9,ip) of the DB to 
be triggered by the fluorescence detector. Q V (E) is the 
flux of high-energy neutrinos depending on the model of 
the extragalactic source of high-energy cosmic rays. E, 
A, T and £1 are the energy of the incident neutrino, area 
over the f.o.v. of the detector, time of data acquisition 
and solid angle around the detector, respectively. 

The interaction probability is given approximately by: 

P int (E,6) = (a uN (E))x(N T ( X )) (3) 

where (<t v n(E)) is the average cross section of the 
neutrino- nucleon interaction and (Nt(x)}, the average 
total number of nucleons per squared centimeter at the 
interaction point in the atmosphere. Nt{x) — %Nax{9), 
where Na is the Avogadro's number and x(9) is the at- 
mospheric slant depth at the neutrino-nucleon interac- 
tion point. 

Considering the Earth's curvature, the slant depth can 
be approximately written as: 

X (9) = J p(H = IcosO + ^^)d\ (4) 

where A is the path along the arrival direction from infin- 
ity until the interaction point in the atmosphere, p is the 
atmospheric density, H the vertical height, I is the dis- 
tance between the interaction point and the point toward 
the particle goes through on Earth (the slant height), and 
9, the zenith angle. The atmospheric depth as a function 
of the zenith angle is shown in Fig. 

As the cross section of neutrinos with ultra-high en- 
ergies is unknown, usually one adopts the extrapolation 
of parton distribution functions and Standard Model pa- 
rameters far beyond the reach of experimental data. In 
this way, one can estimate a value for the cross section 
of the neutrino-nucleon interaction of about 10~ 32 cm 2 , 
for energies around 1 EeV. Some authors say that this 
extrapolati on g ives a neutrino-nucleon cross section that 
is too high |23 but others use models that increase this 
same cross section to typical hadronic cross section val- 
ues [24j. In this work we use the extrapolation of the 
Standard Model cross section. 
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FIG. 2: Number of charged particles as a function of atmospheric depth in g/cm 2 simulated for incident angles of 45°, 55°, 
65° and 75° from left to right and up to down. Inferior lines represent the first and second EAS, and the upper line is the sum 
when the two EAS are superimposed. The energy of the primary neutrino is 0.5 EeV. 



The trigger probability is given by: 

P trig (E,r,6,<p) = rxP had xP db (E,r,6,<p)xZ{E,r) (5) 

where T is the fraction of the time the fluorescence de- 
tector will work (T = 0.1 because the fluorescence detec- 
tor can only operate in clear moonless nights), Phad is 
the hadronic branching ratio of tau decay {Phad = 0.64), 
Pdb(E, r, 6, ip) is the probability of the DB to be seen by 
the detector and r) is the efficiency of the detector. 

We define the probability of the DB to be seen by the 
fluorescence detector as: 

P db (E,r,6, V ) = ^±p- ,tfuj(r,9,<p)<L(E) 

= 1, otherwise (6) 

where, as we have seen in Section^ L(E) is the distance 
traveled by the tau in laboratory frame and u(r, 9, (p), as 
can be seen in the Fig. ^ is the size of the shower axis in- 
side the f.o.v. of the detector. We made a rough approxi- 
mation to account only for the showers moving away from 
the detector, so that J ui(r,6,<p)dip ~ Tru(r,0). Then we 
measure u)(r,6) where the vertical plane containing the 
shower axis passes through the center of the fluorescence 
detector. 

The efficiency Y,{E, r) is the convolution of the energy 
efficiency Y,'(E) and the efficiency depending on distance 



S"(r). As the fluorescence detector has no measurement 
of Yi'(E) yet, but only expectations based on simula- 
tions, we constructed one function based on simulations 
that grows logarithmically from to 1 in the interval 
of neutrino incident energy 0.3 < E (EeV) < 88 [29| . 
£'(£■) = if the energy is lower than the energies in this 
range and H'{E) = 1 if the energy is higher than that. 
S"(r) is a Gaussian distribution centered at the point 
15 km far from the detector, but for simplicity in the 
calculation we used a step function 20 km long centered 
at the same point. It does not make a great difference in 
the final result. 



III. RESULTS AND DISCUSSION 

Because the Auger Observatory Fluorescence Detector 
will be constituted of four detectors with a site seeing of 
180°, we have to multiply the result of Eq. Qby 2. Then 
using Eq. ^ we calculated the number of events which 
can be seen in Tables IIIII and IIVI for different models of 
ultra-high energy cosmic ray flux and in different energy 
intervals. 

From Table IIIII one can learn which is the energy 
interval which is relevant to detect DB events in the 
Auger fluorescence detector. It is, approximately, 0.63 < 
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FIG. 4: Atmospheric depth as a function of the zenith angle. 
The inferior line represents the interaction point depth for 
I = 30 km. 



E (EeV) < 2. The neutrino flux used here for the model 
AGN-95J is an approximation that takes roughly the 
average between the models A and B in reference [32$ . 
Table IIVI shows the expected number of DB events in 



TABLE III: Number of events in the fluorescence detector 
during a period of three years, calculated in different regions 
of the energy spectrum and for different models of neutrino 
flux. TD-92 stands for the model in reference [3(1; TD-96 for 
the reference [H; AGN-95J for H2; WB for ^ and MPR 
for 0. 



Models (p) 


Nx a 


N 2 b 


A 3 C 


N A d 


TD-92 (0) 


260 


450 


180 


70 


TD-92 (0.5) 


2.8 


5.7 


2.8 


1.5 


MPR 


2.4 


3.6 


1.1 


0.3 


TD-92 (1) 


0.14 


0.37 


0.25 


0.2 


AGN - 95J 


0.15 


0.22 


0.07 


0.02 


TD-92 (1.5) 


0.06 


0.15 


0.1 


0.08 


WB 


0.05 


0.08 


0.02 


0.01 


TD-96 (1) 


2.9 x 10" 5 


3.1 x 10" 5 


6.1 x 10" 6 


1.0 x 10" 



"0.20 < £(EeV ) < 0.63 
6 0.63 < £(EeV ) < 2.00 
c 2.00 < E(EeV ) < 6.30 
d 6.30 < E(EeV ) < 200 



three years according to the ultra-high energy cosmic rays 
source models, when a Standard Model extrapolation for 
the ovat is used. 

Based on the simulations shown in Section III Al one 
can have an idea of the longitudinal development of the 
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TABLE IV: Tau neutrino fluxes at the Earth according to 
the corresponding total number of events in a period of three 
years. 



Models (p) 


(GeV" 


1 m" 


2 s _1 sr _1 ) 


Agents (Syrs)" 1 


TD-92 (0) 


2.31 x 10" 


9 T 1 


■ 1 77 


9.6 x 10 2 


TD-92 (0.5) 


5.64 x 10" 


e E- 


1.58 


1.3 x 10 1 


MPR 


2.15 x 10" 


-2 E - 


2 


7.4 x 10° 


TD-92 (1) 


1.07 x 10" 


-9 E - 


1.3 


9.6 x 10 _1 


AGN - 95J 


1.34 x 10" 


-3 E - 


2 


4.6 x 10 _1 


TD-92 (1.5) 


4.24 x 10" 


-io E 


-1.3 


3.8 x 10 _1 


WB 


4.83 x 10" 


4 E- 


2 


1.7 x 10" 1 


TD-96 (1) 


7.51 x 10" 


4 E- 


■2.4 


6.7 x 10~ 5 



DB as a function of the incident angle and energy of 
the primary neutrino. It is the convolution of the terms 
Pdb(E ', r , 9 , tp) and T,(E,r) in Eq. [31 that describes the 
restriction in the energy interval observed in the simula- 
tions and Table ITTT1 For relatively low energies the effi- 
ciency T,(E, r) of the detector is also low and, as can be 
seen in Eq. for energies higher than 2 EeV the factor 
Pdb{E, r, 6>, ip) is too small. 

There is also a restriction for the neutrino incident 
zenith angle. For zenith angles smaller than 55°, be- 
cause the lower atmospheric density, the development of 
the shower is slower and there is a superposition of the 
two EAS. The DB still can be triggered but in a more 
sophisticated way by comparing the longitudinal develop- 
ment of the EAS. Fig.[!)]shows events with incident angles 
equal to 45°. In the simulations on the top of Fig. [S] we 
can observe the two overlapped showers such that they 
look like a single ordinary EAS. Even in that case one can 
see that it is possible to identify a DB because both the 
maximum of the charged particles number and the lon- 
gitudinal development of the shower are not compatible 
with a "single-bang" of the same energy, which is shown 
in the lowest plot of Fig. [5] For zenith angles larger then 
75°, since h (see Fig.^l is too low, it is difficult to detect 
the first EAS. 

We simulated sixty DB events for each different en- 
ergy and angle that we present. The results have a root 
mean square deviation for the maximum of the showers 
of approximately 30 g/cm 2 . As the maximum of the DB 
(top and central plots) in Fig. [S] is around 700 g/cm 2 
and the maximum of the ordinary EAS (lowest plot) is 
around 500 g/cm 2 , it is easy to differentiate between the 
depth of the maximum of a DB and an ordinary EAS in 
the simulations we made. If a higher energy proton will 
masquerade a lower energy DB depends on the accuracy 
of the detector to determine the energy of the primary 
particle. Because of the mean life time of the tau, only 
relatively low energy DB will be superimposed looking 
like an ordinary EAS. So if you detect an ordinary EAS 
profile of relatively high energy (E > 10 EeV), that can- 
not be considered a DB event. 

To study the relation between the cross section and 



the potential of the Pierre Auger Fluorescence Detector 
to detect DB events we wrote the neutrino cross section 
as <7„at = A x (Tsm, where A is a free parameter which 
depends on the model and asM ~ 7.8 x 10~ 36 cm 2 X E 363 
is the Standard Model extrapolation for a v jq which have 
10% accuracy within the energy range 10~ 2 < E(EeV) < 
10 3 when compared with the results of the CTEQ4-DIS 
parton distributions |35| . 

To account for models that gives a cross section differ- 
ent from the extrapolation of the Standard Model cross 
section we just vary A. It is a naive approximation but 
gives an idea of how the number of events can increase 
or decrease if the cross section is not standard. For in- 
stance, if the cross section is of order 10 times higher then 
the standard in the energy range 0.20 < E'(EeV) < 200, 
then even the WB upper bound for the ultra-high energy 
cosmic ray flux predict a rate of more than 1 event in 3 
years in the Pierre Auger Fluorescence Detector. 

The horizontal lines in Fig.[S]are number of events cho- 
sen for comparison (1 event in 3 years is considered the 
lower limit for the rate of events in the Auger Observa- 
tory). The vertical line shows the extrapolated Standard 
Model cross section. The other lines show the number 
of events for different models of ultra-high energy cosmic 
ray flux depending on A. For instance, we can see from 
this figure that assuming the Standard Model extrapo- 
lated cross section, if the Observatory measure a number 
of events around one in three years, models like TD-96(1), 
TD-92(0) and TD-92(0.5) will be excluded, and the more 
compatible with the measurements will be TD-92 (1). 



IV. CONCLUSION 

Taken into consideration neutrino oscillations, one ex- 
pect that one third of the high-energy neutrino flux ar- 
riving at the Earth should be composed of tau neutrinos. 
These neutrinos can interact in the Earth's atmosphere 
generating a double shower event named Double-Bang 
Phenomenon. Many recent works have studied the po- 
tential of the Pierre Auger Observatory to detect horizon- 
tal air showers generated by ultra-high energy neutrinos 
with the surface detector. Here we specifically investi- 
gate the potential of the fluorescence detector of that 
observatory to observe DB events. 

DB events have very particular characteristics in the 
Auger Observatory. Different from the neutrino events 
in the surface detector, DB events do not need to come 
from the very near- horizontal angles. Despite the low 
probability of interaction in the atmosphere, we can also 
have tau neutrinos creating DB events with incident an- 
gles from approximately 55° to 75°. DB events also have 
a lower energy, around 1 EeV, different from the ener- 
gies around 50 EeV and beyond expected for an ordinary 
EAS generated by the highest energy cosmic rays. In 
the range of energy approximately between 0.6 EeV and 
2 EeV a considerable part of the two EAS that charac- 
terize a DB can be detected by the fluorescence detector 
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FIG. 5: Number of charged particles as a function of atmospheric depth in g/cm 2 simulated for an incident angle of 45° 
and depth of the first interaction equal to 24.3km ~ 30g/cm 2 . The first two graphics from the top were simulated using the 
CORSIKA program versions 5.62 (which is valid for zenith angles smaller than 60°) and 6.00 (which is valid even for zenith 
angles larger than 60°) respectively, with the energy summed of the two EAS equal to 0.5 EeV. Inferior lines represent the first 
and second EAS, and the upper line is the sum when the two EAS are superimposed. The graphic in the bottom is an ordinary 
shower simulated with the CORSIKA version 6.00 and energy of 0.5 EeV. 
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FIG. 6: Number of events according to the parameter A of 
the cross section for various models. A = 1.00 corresponds to 
the extrapolation of the Standard Model cross section. See 
text for details. 



and then we have a DB trigger. 

The number of DB events depends on many parame- 
ters like arrival flux, cross section, energy and incident 
direction of the neutrinos, and efficiency of the detec- 
tor. Models like MPR, TD-92(1) and AGN-95J generate 
a number of events of around one in three years in the 
Pierre Auger Fluorescence Detector. This is because the 
energy range where the DB can be detected is very strict. 
For energies less then 0.6 EeV the efficiency of the detec- 
tor is too low and for energies greater then 2 EeV the 
two EAS are too separated. 

Concerning the background, the probability for a pro- 
ton to generate a DB and masquerade the DB generated 
by a neutrino depends on two possibilities: 1) that the 
primary proton interaction generates some fragment that 
will give rise to a secondary shower deep in the atmo- 
sphere with energy higher then the first. 2) that another 
shower created by some independent particle interacts 
deep in the atmosphere masquerading the second EAS of 
the DB. 

In the possibility 1, the second EAS will be created 
by the decay or interaction of the fragment deep in the 
atmosphere. We consider that the primary proton looses 
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roughly half of its energy to the secondary particles that 
constitute the EAS, and so it is very hard that the second 
EAS has more energy than the first one. Now, consider- 
ing that for this high-energies we have a cosmic ray flux of 
the order of 1 particle per km 2 per year and that the only 
particles that could interact deeply in the atmosphere are 
neutrinos, generating the second independent EAS near 
the detector, the chance that the proton and this second 
independent neutrino come from the same solid angle di- 
rection interacting in a time interval of the tau decay in 
the laboratory frame of ft « 131 X M s 1S approxi- 

mately (at most) 1 in 10 s , what exclude the possibility 2. 
The direction of the two EAS can be identified specially 
if two fluorescence detectors trigger the same DB event 
(with only one detector, it must be difficult to know the 
direction of the EAS in the plane that contains the EAS 
and the detector). 

Based on this assumptions, EijE\ can be a good pa- 
rameter to identify DB events if the error and the aver- 
age in the energy measure is within certain specific value. 
The error in the energy measured by the fluorescence de- 
tector will depend mainly on the atmospheric conditions 
but hardly will exceed 50%. For a DB event the situation 
is optimistic because the most important is the relation 
between the energies of the two EAS and this error is 
smaller than the error of the absolute energy of an ordi- 
nary EAS. We can make a conservative estimation of the 
error in the average ratio E2 / E\ considering the error in 
the absolute energy of 50%. This will give a relative error 
to the energy ratio of 70%. So, since EijE\ w 2.67 in 



average as deduced in Section [H] then considering such 
an error we find 95% of the events such that the energy 
ratio EijE\ > 1. Then E^/Ei > 1 is a good parameter 
to identify DB events. 

Despite the fact the DB Phenomenon can be very rare, 
it is very important to be prepared for its detection, spe- 
cially in case the Pierre Auger ground array detect near- 
horizontal air showers which can indicate a sign for elec- 
tron and/or muon neutrinos. Consequently oscillations 
imply a considerable number of tau neutrinos too. With 
such a motivation, the Auger Observatory could calibrate 
its trigger to be more sensitive to energies around 1 EeV 
or increase the number of fluorescence detectors. 

The potential of the DB Phenomenon to acquire valu- 
able information both in particle and astrophysics is ir- 
refutable. For instance, the cross section and flux of the 
ultra-high energy neutrinos are speculative and can be 
investigated with DB events. 
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